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Finite-density quantum chromodynamics (QCD) 


QCD 
Fundamental theory for quarks and gluons 


Neutron star 


e Cold and dense nuclear matter 


@ 2m, neutron star (2010) 


Neutron star merger 
e Gravitational-wave observations (2016~) (image from NASA) 


Path-integral expression of finite-density QCD: 


Zaco(T. n) = [ DA Det(P(A, pa) +m) exp C Svu(). 


quark gluon 


Sign problem: Det(D(A, uqa) 4- m) Z 0 at uq z 0. 
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Sign problem & Complexification of variables 


Consider the path integral: 


ja J «eec sta 


e S|[x] is real = No sign problem. Monte Carlo works. 
e S|[x] is complex = Sign problem appears! 


If S[r] € C, eom S'[x| = 0 may have no real solutions z(t) € R. 


Idea: Complexify x(t) € C! 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) = EL (= + 2 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on A 
(a — 1) 


Yuya Tanizaki (RBRC) Review on Lefschetz-thimble method Aug 2, 2016 @ Plymouth 4/31 


Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (c = 1,2) for the Airy 


integral: 
Ai(a) Son f 2 p am 
= a — —+az). 
= J, 2T p 3 


ng: intersection number of the steepest ascent contour Ke and R. 


Yuya Tanizaki (RBRC) Review on Lefschetz-thimble method 


Gradient flow 


Problem in the multi-dimension 
Im(S) = const. gives (2n-1)-dim. manifolds, instead of n-dim. ones. 


This defines the steepest descent directions, since 


e-XK I 


The flow lines satisfies Im(;S) = const. [Witten, arXiv:1001.2933, 1009.6032] 


Gradient flow Consider 


> 0. 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles Js: (classical eom S’(z,) = 0) 


f d"z e 3) — S (Ke; 3 d'zg 9i 


o 


Jz are called Lefschetz thimbles, and Im[S] is constant on it: 


M =}, dzi(t) _ eo 


1—-—oo dt pe 


= [40 


(Ko, R): intersection numbers of duals K, and R” 


(Ko = {2(0)|2(00) = 253). 


[Pham, ‘83, etc., Witten, arXiv:1001.2933, 1009.6032] 
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MC on one-thimble 


Monte Carlo simulation on one Lefschetz thimble 


Most of the works before LATTICE 2015 are devoted to MC method 
with one-thimble ansatz. 


z= | dz eS) = z= d^z e80, 
" Jo 


[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 
Motivation 


e Within the mean-field approx, this seems to be justified for 
bosonic theories. 


e |t was not known how to take the summation over thimbles. 


It is successful for several models, and a lot of numerical techniques 
are developed. 
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MC on one-thimble 


Relativistic Bose gas: 


S= f at [Jao]? (v? — iol + uea XI 


(Cristoforetti et al., PRD 88 (2013) 051501; Fujii et al., JHEP 1310 (2013) 147; 
Cristoforetti et al., PRD 89 (2014) 114505; Alexandru et al. 1606.02742) 


3L 


|» (number density) 
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1 14 12 13 
u (chemical potential in lattice unit) 


Fujii et al. (JHEP 1310) 
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MC on one-thimble 


Chiral Random Matrix model 


ich(u)X + sh(u)Y ) NG 


S = Nu(X! X +Y'Y) — In det ( jch(u)Xt 4 sh(u)Yt 


CRMT with 1-thimble ansatz with Ny = 2, u/V N = 2. (Di Renzo, Eruzzi, PRD(2015)) 


(cf. Naive CL gives the phase-quenched result. (Mollgaard, Splittorff, 1309.4335) 


Some gauge cooling extends applicability of CL until u/WN € 3 (Nagata et al. 1604.07717)) 
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fermion model 


(0+1)-dimensional fermion model 


List 

e Tanizaki, Hidaka, Hayata, 1509.07146 

e Fujii, Kamata, Kikukawa, 1509.08176, 1509.09141 

e Alexandru, Basar, Bedaque, 1510.03258 

e Alexandru, Basar, Bedaque, Ridgway, Warrington, 1512.08764 
Related studies 


@ 2-dim Hubbard on 1-thimble (Mukherjee, Cristoforetti, 
1403.5680) 


e 0-dim models (Tanizaki, 1412.1891, Kanazawa, Tanizaki, 
1412.2802) 


e Ch. RMT on 1-thimble (Eruzzi, Di Renzo, 1507.03858) 
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fermion model 


One-site Fermi Hubbard model 
One-site Hubbard model: 


Ê = Ufuf, — p(y + fry). 
Fock state gives the number density immediately: 


139 2(e8# + oPQu-U)) 
~ Bon D = 14-29 x PCU)" 


In the zero-temperature limit, 


(YT, Hidaka, Hayata, 1509.07146) 
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fermion model 


Path integral for one-site model 


Effective Lagrangian of the one-site Hubbard model: 
v? 
L= oU + v* [8. — (U/2 + ip + u)] v. 
The path-integral expression is 


B ik d 2 052 
Z24Í — | ay (14 ef Gem)? 78/20 
2nU R eee 


Fermion Det 


Integrand has complex phases causing the sign problem. 
ọ is an auxiliary field for the number density: 
(ñ) = Im(y)/U. 
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fermion model 


Sign problem and Gradient flows at u/U < —0.5 


U f 
Det lo — (» + F + ie) = (1 Ta e AD 


~R. 


Figure: Flow at u/ 


(YT, Hidaka, Hayata, 1509.07146) 
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fermion model 


Flows at —0.5 < u/U < 1.5 


Re(z) 


Figure: Flow at u/U = 0 


Im(z,)/U œ u/U + 1/2. 


Complex saddle points lie on nymf 
(YT, Hidaka, Hayata, 1509.07146) 
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fermion model 


Complex classical solutions 


Classical solutions: 


u[U +1/2 


U 
2m ted B+ > + 2rmT. 


n) CE 


Re (Ss — So) œ aum m?, 
Im S,, ~ 27m (£ + i). 


Classically, Zetassical = X e m. 
m 
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fermion model 


Numerical results 


Results for GU = 30: (1, 3, 5-thimble approx.: Jo, Jo U Ji, and Jo U J41 U Ja ) 
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Necessary number of Lefschetz thimbles ~ GU/(27). 
(YT, Hidaka, Hayata,1509.07146) 
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Relation with complex Langevin method 


e Aarts, 1308.4811, Aarts, Bongiovanni, Seiler, Sexty, 1407.2090 
e Tsutsui, Doi, 1508.04231 

e Fukushima, Tanizaki, 1507.07351 

e Hayata, Hidaka, Tanizaki, 1511.02437 

e Abe, Fukushima, 1607.05436 
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Complex Langevin method 


Complex Langevin has been regarded as a sign-problem solver via 
stochastic quantization (Klauder, PRA 29, 2036 (1984), Parisi, PLB 131, 393 (1983)): 


dz(0) sS 
0: Stochastic time, 7: Random force (n(0)n(0')),, = 20(0 — 0"). 
Ito calculus shows that 


x (Olen(0)) = A(O" (z,(8))), — (0 (24(8)) S (27(0)))n: 


If the I.h.s becomes zero as 0 — oo, this is the Dyson-Schwinger eq. 
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Relation between CL and LT? 


Both methods relies on complexification, but not much is known for 
their relations. 

CL and LT looks similar, but they are still different: 

(U(1) link model , S = —8cos(z) — In[1 + &cos(z — iu)]) 


x 


k= 1/2 
(Aarts, Bongiovanni, Seiler, Sexty, 1407.2090) 
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Relation with complex Langevin 


Semiclassical incorrectness of CL method 


If A is small enough, we can show a sufficient condition for incorrect 
behaviors of CL method. 
(Hayata, Hidaka, YT, 1511.02437) 


Assume that CL method is correct, then 


(Oleh = 7 DKR") f dze sola), 


Since A < 1, 


SM NECEM » G0): 
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Relation with complex Langevin 


Semiclassical inconsistency 


In the semiclassical analysis, one now obtains (for dominant saddle 
points) 


= (Ko, R”) Qh e SG)/h. 
Z S" (zo) 


The right hand side can be complex, which contradicts with c, 0! 
(Hayata, Hidaka, YT, 1511.02437) 
We show that the complex Langevin is wrong if 
e There exist several dominantly contributing saddle points, and 
e Those saddle points have different complex phases. 
Open question Connection with the histogram test on CL method? 
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Relation with complex Langevin 


Refine Complex Langevin via thimbles (1) 


Deform the theory so that only one thimble contributes, and apply CL 
(Tsutsui, Doi, 1508.04231) 


Z= f| fe Sax 1 = fu@ + ig(x))e 9) dz. 
One can compute VEV of the original theory using the new one as 


(QU otiginat = Re(O)new ax haven p tQ), 
EP aigue 


(G) quench/ (f) quench IS common for any observables. 
Compute (O)new using CL with "appropriate" g. 


Open question What g should be chosen? 


Yuya Tanizaki (RBRC) Review on Lefschetz-thimble method Aug 2, 2016 € Plymouth 23 / 31 


Relation with complex Langevin 


Refine Complex Langevin via thimbles (2) 


Perform the reweighting by attaching complex phases of thimbles to 
CL distribution (Hayata, Hidaka, YT, 1511.02437) 
Test on one-site fermion model 


BU =30 &4/U —0 
* --- exact 
HEH w/o phase 


1H w/ phase 


: J -05 00 05 To 15 
Š 2 u[U 


=§ 


0 
BRe(z) 
Clear improvement, but there's unknown systematic error. 
Open question Can we justify and make it rigorous? 


Yuya Tanizaki (RBRC) Review on Lefschetz-thimble method Aug 2, 2016 @ Plymouth 24/31 


Simulation on multiple thimbles 


Simulation on multiple thimbles 


List 
e Alexandru, Basar, Bedaque, Ridgway, Warrington, 1512.08764 
e Alexandru,Basar, Bedaque, Vartak, Warrington, 1605.08040 
Related studies 
e Alexandru, Basar, Bedaque, Ridgway, Warrington, 1604.00956 
e Alexandru, Basar, Bedaque, Ridgway, Warrington, 1606.02742 
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Simulation on multiple thimbles 


Possible concerns for practical applications 


Interference among Lefschetz thimbles is very important for our 
interest (especially when fermion exists). 


This means that we must ... 
e Find all contributing complex saddle points, 
e Construct Lefschetz thimbles around those saddle points, 
e Evaluate integration on each Lefschetz thimbles, and 
e Sum up those results. 
We need some machinery to do them automatically. 
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Simulation on multiple thimbles 


Idea for multiple thimble simulation 


d Os 
Deform the original cycle IR” by the gradient flow, E E (5) 


(Alexandru, Basar, Bedaque, Ridgway, Warrington, JHEP (2016)) 


Im(z) 
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Simulation on multiple thimbles 


Idea for multiple thimble simulation 


À . d oS 
Deform the original cycle IR” by the gradient flow, E E (5) 
Z 
(Alexandru, Basar, Bedaque, Ridgway, Warrington, JHEP (2016)) 
Im(z) 
I(T) 
Re(z) 
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Simulation on multiple thimbles 


Formulation 
Let us fix a flow time T, and define 


JU e far xr) € C” da) = (5): A054) =2 2 | i 


Aa 


dt Oz 


By construction, z(T;-) : IR^ Š J (T) and 


] seo = i dz e-50) 
n I(T) 


i d"z det (22) e- S62). 
n Ox) 


=>One can do usual Monte Carlo + reweighting by regarding 


Seer(t) := S(z(T;2)) — In L (2) 


oxi 


as the effective classical action. 
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Simulation on multiple thimbles 


Real-time dynamics 
This method is applied to Schwinger-Keldysh path integral for 


u EU 
s- fa (s 52 2 


0.3 È a LEN EE Re[Cr (t.t)] 
= EN — Im[Gr (tt) é $ 


0.0 0.5 1.0 1.5 2.0 


Feynman propagators at 3 = 0.8. Tgow = 0.2. (Alexandru,Basar, Bedaque, Vartak, 


Warrington, arXiv:1605.08040) 
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Summary & Conclusion 


Summary and Conclusion 


e Lefschetz-thimble method is helpful to analyze structures of sign 
problems. 

e Many Lefschetz thimbles can contribute. Especially, interference 
among them will play an important role for physical observables. 

e Dynamics in complexified space is complicated. 
Comparison among one-thimble ansatz, complex Langevin, and 
saddle point analysis gives us a good insight. 

e Recent developments may enable us to study nonperturbative 
field theories with the sign problem. 
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